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Abstract. A new method for a simultaneous search for clusters of galaxies in X-ray photon maps and optical 
galaxy maps is described. The merging of X-ray and optical data improves the source identification so that a large 
amount of telescope time for spectroscopic follow-up can be saved. The method appears thus ideally suited for 
the analysis of the recently proposed wide-angle X-ray missions like DUO and ROSITA. As a first application, 
clusters are extracted from the 3rd version of the ROSAT All-Sky Survey and the Early Date Release of the Sloan 
Digital Sky Survey (SDSS). The time-consuming computations are performed within the German Astrophysical 
Virtual Observatory (GAVO). On a test area of 140 square degrees, 75 X-ray clusters are detected down to an 
X-ray flux limit of 3 - 5 x 10~^^ ergs"^ cm"^ in the ROSAT energy band 0.1-2.4 keV. The clusters have redshifts 
z < 0.5. The survey thus fills the gap between traditional large-area X-ray surveys and serendipitous X-ray cluster 
searches based on pointed observations, and has the potential to yield about 4,000 X-ray clusters after completion 
of SDSS. 
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1. Introduction 

Clusters of galaxies trace the peaks of the matter distri- 
bution and their dynamics on cosmic scales (see review in 
Borgani & Guzzo 2001). The corresponding cluster num- 
ber density fluctuations can be measured up to Gigaparsec 
scales where they display a simple Gaussian random field 
with a high statistical significance (Schuecker et al. 2002, 
2003a). Therefore, the related mean spatial abundance, 
the two-point correlation function, and the power spec- 
trum give a useful summary of these fluctuations. They 
are well fit by a flat, low density Cold Dark Matter model 
without any significant contribution of dark energy others 
than the cosmological constant (Schuecker et al. 2003b). 
The data further support a simple scale-invariant biasing 
scheme although the details are still discussed controver- 
sal (e.g., Bahcall 1988, Collins et al. 2000, Schuecker et al. 
2001, Zandivarez, Abadi & Lambas 2001). 

The usuability of rich clusters for cosmological inves- 
tigations is mainly related to their simple statistical prop- 
erties and to their comparatively simple physical struc- 
ture (e.g. Kaiser 1986, Sarazin 1988, Neumann & Arnaud 
2001). This is especially the case for X-ray clusters of 
galaxies where the X-ray flux limit provides a clean mea- 
sure of the cluster selection function. A crucial point is the 
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reliability of cluster mass estimates (Schindler & Miiller 
1993, Evrard, Metzler & Navarro 1996) where evolution- 
ary effects could further complicate the tests (Vikhlinin 
et al. 2002). However, for not too distant systems the ef- 
fects are found to be small (e.g. Gioia et al. 2001, Rosati, 
Borgani & Norman 2001) so that a proper calibration with 
detailed X-ray observations appears to be possible and is 
already in progress (H. Bohringer et al., in preparation). 

Of fundamental importance for every precise cosmo- 
logical test is the construction of large, homogeneously se- 
lected and thus representative samples. Unfortunately, the 
time-consuming spectroscopic follow-up observations of 
the cluster candidates selected from X-ray data hampered 
the construction of X-ray cluster samples with sizes of the 
order of 10^ or more (see review in Edge 2003). A sim- 
ilar problem is expected for the recently proposed large- 
area X-ray missions like DUO and ROSITA. From the 
construction of the Northern ROSAT All-Sky (NORAS, 
Bohringer et al. 2000) and ROSAT ESO Flux-Limited X- 
ray (REFLEX, Bohringer et al. 2003) samples we learned, 
however, that a careful pre-selection of cluster candidates 
can lead to high success rates during the spectroscopic 
follow-up observations and can thus save a large amount 
of telescope time (see also Bohringer et al. 1998, 2001, 
Guzzo et al. 1999). 

In the following section we give a description of a 
general likelihood method for the detection of clusters of 
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galaxies (Sect.|21) which allows a simultaneous analysis of 
multiwavelength data. In Sect. Owe apply the method to 
the X-ray photon distributions as given by the 3rd version 
of the ROSAT All-Sky Survey (RASS-3, Voges et al. 1999). 
A new source detection in RASS is necessary because it 
will turn out that - because we combine X-ray data with 
much deeper optical data - the new method can go signifi- 
cantly deeper compared to the standard analyses of RASS. 
For the measurement of the source extent, hardness ratio 
and other important X-ray quantities the Growth Curve 
Analyses (GCA) method of Bohringer et al. (2000, 2001) 
is used. In Sect.0]we apply the same likelihood method to 
the galaxy distribution as given by the Early Date Release 
of the Sloan Digital Sky Survey (SDSS, York et al. 2000) 
using the r' band where completeness and sample sizes are 
highest down to the faintest magnitudes and which serves 
as a 'reference' passband where most methods using SDSS 
data are first illustrated (future applications will, however, 
include all five SDSS bands). The similarity of the source 
detection in X-rays and cluster search in the optical allows 
a straightforward merging of the two results (Sect.EJ. 

In a test area of 140 square degrees, 75 clusters 
are detected down to the nominal flux-limit of 3 — 5 x 
10-13 erg s"^ cm"2 in the ROSAT energy band 0.1 - 
2.4 keV, that is, 5-10 times deeper than REFLEX. The 
measured parameters of the clusters are described in 
Sect. El Extrapolated to a final SDSS area of 7,000 square 
degrees, our method is thus expected to yield about 4,000 
X-ray clusters. To test the reliability of our method, we 
compare in Sect. [3 our X-ray cluster sample with various 
SDSS and RASS-based catalogues. Further tests with re- 
ahstic simulations are postponed to a second paper. 

It will be seen that a detailed analysis of the RASS and 
SDSS data requires massive computing power. Moreover, 
our method allows the inclusion of further high-resolution 
maps from, e.g., temperature anisotropics of the cosmic 
microwave background radiation, centimeter or millimeter 
(Sunyaev-Zel'dovich) maps etc. The combined analysis of 
such large databases is a typical example where virtual ob- 
servatories can help with a fast and easy data handling and 
a distribution of the time-consuming computations over 
many computer systems (grid computing). The present 
investigation was thus planned already from its beginning 
to be performed within the German Astrophysical Virtual 
Observatory (GAVO). 

All computations assume a pressureless Friedmann-Le- 
maitre world model with the Hubble constant Ho in units 
of h = iJo/(100 kms^^ Mpc^^), the normalized cosmo- 
logical constant JIa — 0.73, and the density parameter 
Vtm = 0.27, but the resulting cluster catalogues do not de- 
pend on the chosen values of the cosmic energy densities. 

2. A general method for cluster detection 

In the following we want to define likelihood functions ap- 
propriate for the detection of clusters in point patterns. 
For the analysis a quantity Lv{xi, . . . ,xn) dxi ■ ■ ■ dx^ 
is defined which gives the probability of finding exactly 



one point in each of the infinitesimal volume elements 
dxi, . . . ,dxN centered a.t Xi, . . . ,xn. We assume that this 
probability is Poissonian with the intensity parameter A. 
The presence of a cluster suggests that A is closely related 
to the projected cluster profile (plus background) and is 
thus a function of x. The probability to find a point in 
a given volume element is thus guided by A (a;), and is 
independent of the probability to find somewhere else an- 
other point. Because of this independence property, Ly is 
simply the product of N terms, each of the form X{xi)dxi. 

However, the test for the presence of a cluster is still 
incomplete because we did not include the fact that no 
further points were observed outside the N infinitesimal 
volume elements. This probability is given by the void 
probability function exp(— J X{x)dx) where the integra- 
tion extends over the volume V reduced by the sum of 
the infinitesimal volume elements, dxi, which can be ne- 
glected due to their small size so that, in total, the inte- 
gration extends over the complete volume V. Therefore, 
an appropriate likelihood function for cluster search is of 
the form 

N 

\iiLv ^ ~ \{x)dx + lnA(a;i), (1) 

which assumes high values when the model distribution 
A(a;) fits the observed pattern of N points. Often, addi- 
tional properties of each point, like its magnitude, colour 
or energy are available. This information can be directly 
incorporated, extending the dimensionality of by a fur- 
ther space which includes the additional marks m, 

XnLvxM ~ / X{x,m)dmdx + > \nX{xi,mi) . 

Jv Jm ^^-^ 

(2) 

A mathematical exact derivation of ^ is given in 
Schuecker & Bohringer (1998) where the correctness of 
this model was verified with a large number of Monte- 
Carlo experiments. Such additional tests are necessary be- 
cause sometimes likelihood estimates can yield biased re- 
sults. Different versions of this matched filter can be found 
in e.g. Postman et al. (1996), Kepner et al. (1999), Kim 
et al. (2002). The corresponding detection probability is 

P,£_iS.l_exp(-l„A), (3) 

with the likelihood ratio L/Lq giving a measure of the 
probability for the presence of a cluster normalized to the 
probability for the absence of a cluster (blank field) . In the 
following sections we will specialize Eq. to the analysis 
of the spatial distribution of RASS photons, and Eq. Q to 
the analysis of the spatial and magnitude distributions of 
SDSS galaxies. The likelihood for a cluster present in both 
X-rays and optical is then simply given by the point-wise 
product of the two likelihood maps. 
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3. Cluster detection in RASS-3 

For the cluster detection in RASS-3 we assume A to be 
proportional to the /ocaZ countrate c(0), giving the numbei 
of X-ray photon events observed in a unit time interval in a 
specific energy band per arcsec^ at the angular coordinate 
9 measured (in arcsec) from the centre of the region to be 
tested for the presence of a cluster. If t{0) is the exposure 
time (in seconds) of RASS-3 at then Eq. |^ translates 
into 



InL 



N 



£ec{e)t(e) + ^ in [c(6>i)t(6>i)] 



(4) 



For the local countrate we assume the model c{0) = c + 
KP{0/9c): a local background countrate c plus a cluster 
with the (convolved) apparent cluster profile P{6/9c) and 
the projected core radius 9c- With the normalization 
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and the definition of an effective exposure time 
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we get the likelihood for a cluster X-ray source in units of 
the likelihood for the local background without any source, 
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where Lq is given by Q for A = 0. From {Tj) and the 
maximization condition d\n{L/Lo)/dA = the following 
two equations are obtained, 
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(8) 



(9) 



which are used for cluster detection and countrate de- 
termination in the following manner: For a chosen Right 
Ascension, Declination and redshift, the cluster profile in 
form of a standard /3 profile (Cavalicre & Fusco-Femiano 
1976) with the core radius Tc and the slope parameter 
(3 is transformed from metric to angular scales 9c and 
convolved with the average RASS-3 pointspread function 
of the X-ray energy band to give P{6/9c)- The expo- 
sure times t{9) are computed for each position covered by 
RASS-3 and yield the effective exposure time from Eq. 
the countrate from Q, and thus the likelihood ratio from 
(PJ). The filter is shifted over the three-dimensional survey 
volume so that sky coordinates, redshifts and countrates 
which simultaneously maximize Eq. (jSJ indicate the pres- 
ence of an X-ray cluster with the detection probability 
(Eq.EJ. 
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Fig. 1. X-ray /3 profile of galaxy clusters with the slope pa- 
rameter /3 = 2/3, the core radius 150 kpc at different red- 
shifts, convolved with the point spread functions of RASS-3 
for the energy 1.25 keV and normalized according to Eq.|3 The 
RASS-3 point spread function closely resembles the curve la- 
beled z = 1.00. 

Distant and/or compact X-ray clusters might appear 
point-like in RASS-3. For such sources the countrates con- 
verge to quite stable values when the cluster filter slides 
towards higher redshifts. The same is true for real point- 
like X-ray sources. Therefore, the cluster-finder described 
above can also be used for the detection of point sources 
which can be identified by their high value of the esti- 
mated redshift (usually the maximum z of the tested red- 
shift range). 

Another important quantity is the number of source 
counts, i.e., the number of X-ray photons above the X-ray 
background: 



Aph ^ At. 



off ■ 



(10) 



The number of source photons gives a further handle to 
evaluate the significance of the detected X-ray source and 
the expected errors of X-ray countrate and fiuxes plus an 
additional error of the background. 

3.1. Application to RASS-3 

For the cluster detection on RASS-3 we choose a test area 
of 140 square degrees with 0.0 < R.A. < 56 deg and 
-1.25 < DEC. < -fl.25 deg (part of the SDSS Early 
Data Release). Our analysis is restricted to the ROSAT 
energy range 0.5-2.0 keV where the contrast between X- 
ray background and cluster is highest (Bohringer et al. 
2000). GAVO provides one photon event file and one expo- 
sure map of the complete sky based on RASS-3. In the test 
area and the selected energy band we extracted 180,466 
X-ray photons. The median exposure time is 369 seconds 
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Fig. 2. Maximum likelihood contours (la, 2a ...) for RASS-3 in the energy range 0.5-2.0 keV. Crosses mark objects found 
in the merged ROSAT Bright and Faint (RBF) Source Catalogues. Bright X-ray sources at (R.A. = 20.7°, DEC = -0.70°), 
(R.A. = 33.65°, DEC = -0.775°), (R.A. = 40.675°, DEC = -0.025°) and (R.A. = 54.20°, DEC = 0.60°) have been re- 
moved from the RASS-3 photon maps. Note that RBF sources are also shown slightly outside of the SDSS Declination range 
[-1.25,-1-1.25] degrees. 



which is close to the median exposure time of 388 seconds 
of the total RASS-3. We assume a standard /3 profile with 
a fixed shape parameter of /3 = 2/3 and a fixed core ra- 
dius of Tc = 150/i~^kpc. The profiles are convolved with 
the spatially averaged point spread function of RASS-3 
at 1.25 keV (Boese 2000). A sample of normalized cluster 
X-ray profiles (see Eq.O is plotted in Fig.^ 

For the matched filter only the photons located in a 
quadratic detection cell with a total area of 1.0 square de- 
gree are used (typically 1,200 photons in the energy band 
0.5-2.0 keV). This detection cell is shifted by increments 
of 1.5 arcmin along R.A. and DEC. over the test area. 
The size of the increment is comparable to the half power 
radius (96 arcsec) of the point spread function of RASS-3 
at 1 keV. In redshift direction the increment is Az = 0.02 
in the range 0.02 < z < 1.0, consistent with the error 
estimates obtained with the optical data (Sect.[5j. Future 
investigations will have more computing power and will 
decrease the increment sizes in all three directions signifi- 
cantly. The local background countrate is estimated with 
formal Poisson errors < 3% by the median of the exposure- 
weighted countrates obtained in a large number of subcells 
per detection cell. Small background errors are of funda- 
mental importance for the detection of very faint X-ray 
sources. 



The comparatively large size of the detection cell 
yields stable background estimates. However, bright point 
sources with countrates of more than about 0.5 s~^ affect 
the source detection and background determination sig- 
nificantly. They are detected in a first run of the matched 
filter and then removed from the RASS-3 photon event 
file. In the present case we deleted the photons of on av- 
erage one bright point source per 40 square degrees. After 
their removal, an almost unbiased search for faint sources 
can be performed in a second run (see Figs. l2l4ll . 

Figure|5] compares the likelihood contours as obtained 
with the matched filter and the source positions as pub- 
lished in the combined ROSAT Bright and Faint (RBF, 
Voges et al. 1999) Source Catalogues (crosses). The red- 
shift direction is supressed by replacing the likelihood val- 
ues obtained along the z direction by their maximum 
ln(L/Lo) value. The latter likelihood together with the 
R.A. and DEC. are shown. Not shown is the redshift where 
the maximum likelihood is obtained and the correspond- 
ing countrate. In the following we will call these two- 
dimensional ln(L/Lo) arrays maximum likelihood maps. 
The maximum likelihood contours start at ln(i/Lo) = 1.0 
which corresponds to a detection probability of P = 0.63, 
i.e., at the ~ Icr detection limit of RASS-3. 

A detailed comparison of the results in Sect .[T] shows 
that the matched filter finds more X-ray clusters than pub- 
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Fig. 3. Maximum likelihood contours based on RASS-3 X-ray photons (1st and 4th panels, la, 2a ... contours), SDSS galaxies 
(2nd and 5th panels, > 10a) and the combined maximum likelihood contours of RASS-3 and SDSS data (3rd and 6th panels, 
> lOcr). Crosses mark the position of the deepest X-ray cluster samples available sofar (REFLEX-2, X-ray flux limit 1.8 x 
lO"^'^ ergs~^ cm~^). Squares mark the position of the X-ray clusters of the final sample. The structure at R.A. = 30deg and 
DEC = -1-0.5 deg has a too low likelihood in the SDSS galaxy distribution and is thus not included in the final catalogue. 
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Fig. 4. Maximum likelihood contours. Symbols as in Fig.^l Squares mark the position of the X-ray clusters of the final sample. 



lished in RBF. This is due to the fact that the matched 
filter uses both the RASS and the SDSS data which helps 
to go below the usual detection threshold of the RBF. 
We also have a preliminary MPE-internal list of RASS- 
3 sources obtained with the same detection methods as 
used for the extraction of the RBF sources on RASS-2, 
but without a detailed visual screening of the detections. 
The quantitative comparison in Sect. includes this pre- 
liminary source list. 



Figuresini and 0] show (again) the maximum likelihood 
contours obtained with the RASS-3 survey in the test area, 
but now in comparison to maximum likelihood maps ob- 
tained with the SDSS data and maximum likelihood maps 
obtained from the combination of RASS and SDSS data. 
The determination of the SDSS maximum likelihood maps 
is described in Sect.01 the combination of RASS and SDSS 
maximum likelihood maps and preparation of the final 
cluster sample in Sect.^l REFLEX-2 clusters and the X- 
ray clusters of our final RASS/SDSS sample are shown as 
crosses and squares, respectively, on the combined maps. 
The REFLEX-2 clusters are clearly visible as the most 
significant clusters and illustrate that the present method 
reaches much deeper X-ray flux limits compared to exist- 
ing RASS-based cluster catalogues. 



4. Cluster detection in SDSS 

The detection of clusters in SDSS is based on Eq. jSJ. Our 
analytic treatment follows Sect. Eland is described in detail 
in Schuecker & Bohringer (1998). The final equations are 
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(11) 



,(12) 



which correspond to the Eqs. H8I9|I derived for the X-ray 
data. In Eqs. (|11II12I) . A is now the optical richness which 
is corrected as described in Schuecker & Bohringer (1998), 
h{m) the average number counts of 'background' galax- 
ies with apparent magnitude m, P(-) a King-like profile, 
and (/)(•) the apparent Schechter luminosity function with 
characteristic apparent magnitude m*. Compared to the 
optical richness, we regard the X-ray luminosity as more 
directly related to the total gravitating cluster mass (see, 
e.g., Borgani & Guzzo 2001). Therefore, in the following 
we do not much concentrate on a discussion of the optical 
richnesses. 

4.1. Application to SDSS 

For the optical cluster search, SDSS galaxies with 
extinction-corrected model magnitudes r' < 21.0 mag are 
used. Above this magnitude limit a robust star/galaxy 
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separation is expected (Scranton et al. 2003). In our test 
area we have 2,086,654 SDSS galaxies. A detailed descrip- 
tion of our galaxy selection, especially the choice of SDSS 
catalogue flags, is given in Sect. 4 of Popesso et al. (2003). 

For the optical cluster search a truncated King pro- 
file with a core radius of 150/i~^kpc and a cutoff radius 
of 1 /i-i Mpc is used. The magnitude distribution of the 
background galaxies is obtained in detection cells with an 
area of 5.0 x 2.5 square degrees (about 180,000 galaxies 
per detection cell) . Although this area might be regarded 
as comparatively large for the determination of the local 
galaxy background, we prefer to work with a larger de- 
tection area because the RASS-3 clusters are mainly at 
comparatively small redshifts around z « 0.15 (Sect.[7|l 
where a good measurement of the bright end of the galaxy 
luminosity function is needed. A larger area is also better 
for the detection of nearby less rich systems (Sect.O. 

For the cluster luminosity function we assume a 
Schechter function with the characteristic magnitude 
M*, ~ — 21.43 + 5\ogh and the faint end slope of a = 
—0.90 as obtained from the composite luminosity func- 
tion of SDSS clusters (Goto et al. 2002b). For the con- 
version of absolute into apparent magnitudes we assumed 
the fC-corrections for the early-type galaxies as given in 
Fukugita, Shimasaku & Ichikawa (1995). The richness esti- 
mates are obtained 3.0 mag below the apparent Schechter 
characteristic magnitude. 

In order to allow a direct combination of the likelihood 
contours, the values of the increment parameters in Right 
Ascension, Declination, and redshift used for the analysis 
of the SDSS data are the same as for the X-ray data. The 
maximum likelihood contours are shown in Figs. 13141 

5. Combined cluster detection in RASS-3 and 
SDSS 

The similarity of the likelihood filters applied in X-rays 
and the optical, and the identical sampling of the three- 
dimensional likelihood parameter space (R.A., DEC, z) 
makes the combination of the results obtained with both 
data sets by point-wise multiplication of the likelihood 
values in principle straightforward. However, it turned out 
that the cluster redshift estimates obtained from the SDSS 
data are much better compared to the estimates obtained 
from the RASS data (see Tabs. 2-3). This is mainly caused 
by the limited angular resolution of RASS. Therefore, we 
define the final angular coordinates R.A. and DEC. of a 
cluster by the location of a peak in the multiplied RASS 
and SDSS maximum likelihood map (Figs.|30). The red- 
shift of a cluster is defined by the location of a peak along 
the z direction in the likelihood map of the SDSS data at 
the final (R.A.,DEC.) coordinates. 

The naive multiplication of maximum likelihood maps 
from X-ray and optical data and the subsequent peak anal- 
yses in the combined map can lead to the detection of 
spurious X-ray clusters when a strong signal from one en- 
ergy band is combined with a spurious signal from the 
other band. The resulting systematics can be reduced by 
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Fig. 5. Cumulative X-ray cluster number counts of the 
RASS/SDSS clusters (histograms) for a log- likelihood mini- 
mum of 15 applied to the SDSS data (continuous line), for 
25 (lower dashed line), and for 5 (upper dashed line). The 
RASS/SDSS cluster counts are compared with results ob- 
tained with other surveys (points with Icr error bars: with 
RDCS from Rosati et al. 1998 converted to the ROSAT energy 
band 0.1-2.4 keV assuming a Raymond-Smith spectrum with 
2 = 0, ksT = 5keV, and solar abundance, with REFLEX, 
and REFLEX-2). No corrections for variations of the angu- 
lar survey sensitivity (effective survey area) are applied to the 
RASS/SDSS and REFLEX-2 data (see main text). 



analysing only those peaks in the combined likelihood 
maps where the related peaks in the X-ray and in the 
optical maps exceed certain likelihood thresholds. 

The minimum likelihood threshold for the X-ray data 
is determined by our goal to search for X-ray clusters down 
to the detection limit of RASS. In the present investiga- 
tion we thus use a minimum likelihood of ln(L/Lo) = 3.0 
for RASS-3 which corresponds to a detection probability 
of P = 0.95. The resulting error rate in RASS is compen- 
sated by the high minimum likelihood threshold for SDSS 
of ln(L/Lo) = 15. This value is found to approximately re- 
produce the detection of all known comparatively bright 
REFLEX clusters and limits the number of spurious de- 
tections as shown by the flux/number counts in Fig.|Sl The 
determination of cluster X-ray fluxes and related quanti- 
ties is described in Sect.|H| For this preliminary calibration 
of the likelihood threshold the RDCS sample was not in- 
cluded because the counts have comparatively large sta- 
tistical errors in this range. We show the raw RASS/SDSS 
cluster counts uncorrected for the effective survey area in 
order to illustrate more directly the results of the clus- 
ter search. The final calibration will be performed over a 
much larger sky area including a proper weighting with 
the effective survey area. 
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Fig. 6. Fraction of sources with minimum source counts A'ph 
as obtained with the matched filter for the final X-ray cluster 
sample derive in Sect.|^ 
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Fig. 7. Comparison of countrates obtained with the matched 
filter A and with the growth curve analysis ctr(GCA) for the 
final X-ray cluster sample. Not shown are clusters partially out- 
side the test area (1 cluster) and clusters where GCA could not 
find a clear fiattening of the background-subtracted countrates 
beyond the so-called outer significance radius (6 clusters). 



6. Cluster characteristics 

In our test area, 75 clusters fulfill the criteria derived in 
Sect. El (see Tabs. 2-3, squares in Figs. l3l4f) . Their posi- 
tions, source counts, effective exposure times, countrates 
etc. are fixed by their values at the local maxima in the 
combined maximum likelihood distributions. The X-ray 
source counts A'^ph are then determined with Eq. I|l()|l . The 
distribution function of the number of source counts A'ph 
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Fig. 8. Comparison of countrates obtained from RASS-3 with 
countrates as given in the 2nd ROSAT Source Catalog of 
Pointed Observations with the Position Sensitive Proportional 
Counter. The comparison includes all types of X-ray sources. 
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Fig. 9. Comparison of cluster redshifts estimated from the 
SDSS data with spectroscopic cluster redshifts (see Tabs. 2-3). 

obtained with the final X-ray cluster sample (see Sect.[SJ) 
is shown in Fig.|Sl The cumulative distribution shows that 
about 2/3 of the X-ray clusters have at least 10 X-ray pho- 
tons in the energy range 0.5-2.0 keV above the local back- 
ground yielding formal errors of the countrates of about 
30% (plus negfigible errors of < 3% from the background). 

Fig-El compares the countrates obtained with the 
matched filter A and with GCA for our final cluster sam- 
ple. For countrates A > 0.05/s, no significant differences 
between A and GCA countrates are seen. For A < 0.05 s~^, 
the matched filter gives on average 15% higher countrates 
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Tab. 1. Comparison of countrates A obtained from RASS- 
3 and as given in the 2RXP catalogue. Col. 1: Minimum 
countrate (RASS-3). Col. 2: Number of X-ray sources. Col. 3: 
Average difference (ctr(RASS3) minus ctr(2RXP)). Col. 4: 
Standard deviation of the countrate differences. See also Fig.|H| 
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compared to GCA. We attribute this systematic to the 
fact that the matched fiher integrates the source pho- 
tons out to the edge of the detection cell (formally the 
integration extends to infinity) whereas the GCA inte- 
grates the source photons out to the so-called outer signifi- 
cance radius of a source. Therefore, the former method has 
the danger of slightly overestimating the countrate, while 
the latter method does not count the outer wings of the 
sources and thus slightly underestimates the countrate. 
The differences between these two integration schemes in- 
crease with decreasing countrate as seen in Fig.[71 

GCA further gives for each cluster an additional es- 
timate of the missing flux in order to transform apera- 
ture to total X-ray countrates. For example, at a flux of 
3 X 10^^^ ergs^^ cm^^ the average of the latter correction 
is about 10%. For the matched filter numerical simulations 
are in preparation to measure the appropriate corrections. 

In order to get more information about the relia- 
bility of the matched filter countrates, we also com- 
pare with the countrates obtained from ROSAT PSPC 
pointed observations as published in the 2nd ROSAT 
Source Catalog of Pointed Observations with the Position 
Sensitive Proportional Counter (2RXP, Fig.lHJ, indepen- 
dent whether the target is an X-ray cluster or another X- 
ray source. The comparison is not unrealistic because the 
majority (80%, see below) of the detected faint X-ray clus- 
ters appear point-like in RASS-3a. Averaged systematic 
and random errors of the countrates obtained with RASS- 
3 and ROSAT PSPC pointed observations are shown in 
Table 1. We define a characteristic lower A limit by the 
value where the error of the countrate has the same size 
as the countrate itself. In this case, the comparison sug- 
gests a lower limit of about A = 0.04 s~^. 

Similarily, the systematic errors of the A countrates 
as estimated from the comparison with pointed observa- 
tions range from an underestimation of the countrates of 
(ta/A = 3.9% for A > 0.1 s~^ to an underestimation of 
9.3%atA>0.04s~^. Below this countrate limit the com- 



parison suggests a possible overestimation of the coun- 
trates. One could attribute this error to non-negligible 
contributions by neighbouring X-ray sources. However, 
the systematics are of the order of only 10% compared 
to the random errors at these very low countrates so that 
one should not overestimate the significance of this effect. 

After the countrate/flux conversion of the cluster can- 
didates (see below) we obtain the cumulative flux-number 
count histogram shown as the continuous line in Fig.El 
Only the raw counts uncorrected for the effective survey 
area are shown so that the output of the cluster search 
is seen more directly. The comparison shows that our sur- 
vey fills the gap between traditional large-area surveys like 
REFLEX and faint serendipitous surveys like the RDCS 
(Rosati et al. 1998) which have comparatively large sta- 
tistical errors in this range. A lower limit of the flux limit 
can be obtained from the location where the histogram 
significantly flattens towards faint fluxes (almost indepen- 
dent of the threshold likelihood of the SDSS data). We 
estimated a nominal flux limit of 3 — 5 x 10~^^ erg s~^ cm^ 
(0.1-2.4keV). Below this flux limit a strong incomplete- 
ness is expected caused by, e.g., surface brightness. A more 
quantitative analysis needs a larger test area as provided 
by the SDSS Data Release 2 (coming soon) and a more 
detailed work on simulated data (P. Schuecker et al., in 
preparation) . 

The comparison of the z estimates obtained with the 
optical data with spectroscopic redshifts shows no signif- 
icant systematic deviations and underlines the quality of 
our maximum likelihood cluster redshift estimates (see 
Fig.ini and Table 2). Note that X-ray clusters up to red- 
shifts of at least z = 0.5 can be detected with our method 
with average redshift errors of about = 0.03. 

Tables 2 and 3 give a summary of the X-ray clusters se- 
lected from the combination of RASS-3 and SDSS data. In 
Table 2, column 1 contains a running number of the X-ray 
clusters. Columns 2 and 3 give the cluster sky coordinates 
(in degrees. Equinox 2000.0) computed by weighting with 
the combined maximum likelihood maps of the X-ray and 
optical data. Column 4 gives the column density of neu- 
tral galactic hydrogen in units of 10^" cm^^ as obtained 
from the 21cm observation of Dickey & Lockman (1990) 
and Stark et al. (1992). Column 5 gives the unabsorbed 
X-ray flux in units of 10^^^ erg s~^ cm^^ in the ROSAT 
energy band 0.1-2.4keV determined for an assumed tem- 
perature of 5keV and solar metallicity at redshift zero. For 
the countrate/flux conversion. Tab. 2 in Bohringer et al. 
(2003) was used. Columns 6 and 7 give the maximum like- 
lihood values obtained from RASS-3 and SDSS at the fi- 
nal (R.A.,DEC.) positions, columns 8 and 9 the estimated 
cluster redshifts obtained from the X-ray and optical data, 
respectively. The cluster redshifts obtained from RASS-3 
with values smaller than z = 1 indicate the detection of 
an X-ray extent. A more precise measure of source ex- 
tent is given in Tab. 3 (see below). Column 10 gives the 
spectroscopic cluster redshift as found in the NASA/IPAC 
Extragalactic Database (NED) and column 11 the number 
of cluster galaxies used to determine the spectroscopic red- 
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Fig. 10. RS17 at z = 0.1092 with the X-ray flux S = 1.4 x 10 
at z = 0.2745 with S = 1.3 x 10*^^ erg cm"^ 

shift. The cluster redshifts are obtained from either a pub- 
Hshed cluster redshift or by the median of several galaxies 
with measured redshifts. Galaxies are rejected when they 
fall outside a 3000 km/s velocity interval around the clus- 
ter redshift. In some cases, the remark column in Table 3 
includes additional information on z in form of the cluster 
redshift zq from Goto et al. (2002a). Column 12 indicates 
whether the X-ray cluster is not included in several se- 
lected X-ray and cluster source lists. The notation is as 
follows: 

V2: The ROSAT Bright and Faint Sources Catalogues 
based on RASS-2 (Voges et al. 1999); V3: a preliminary 
(unscreened) MPE internal source list based on RASS- 
3. For the detection of counterparts we assume a search 
radius of 10 arcmin. 

G: The SDSS cluster catalogue published in Goto et 
al. (2002a). For the detection of counterparts we assume 
a search radius of 4 arcmin. 

B: The SDSS cluster catalogue published in Bahcall 
et al. (2003). Note that this list contains a subsam- 
ple of a merged list of SDSS clusters obtained with a 
Hybrid Matched Filter (Kim et al. 2002) and with a color- 
magnitude red-sequence maxBCG technique (Annis et al. 
2003). The Bahcall et al. list has a threshold richness cut 
below the Abell richness class and accepts only clus- 
ters with redshifts 0.05 < z < 0.3. For the detection of 
counterparts we assume a search radius of 4 arcmin. 

Table 3 gives more information about the selected X- 
ray clusters. Columns 1-3 are the same as in Table 2 and 



of X-ray Clusters 




ergs"^ cm"^ in the ROSAT energy range 0.1-2.4 keV. RS28 



contain the cluster number and the sky positions. In col- 
umn 4 the source counts in the ROSAT energy range 0.5- 
2.0 keV at the final sky position are given. Columns 5 and 
6 contain the countrates determined with the matched fil- 
ter and with the CCA method, respectively. The formal 
1(7 errors include the contribution from source count and 
background. In column 8 the difference between the ob- 
served hardness ratio and the theoretically expected hard- 
ness ratio as calculated for a 5 keV cluster and normalized 
to the error of the hardness ratio measurement is given. In 
column 9 the negative logarithmic Kolmogorov-Smirnow 
probabilities for a point source are given. The hardness ra- 
tios and the extent probabilities are calculated with CCA. 
The optical richness can be found in column 9. Special re- 
marks concerning the X-ray clusters are summarized in 
column 10. 

The X-ray/optical overlays in Figs. llOlTTl show typi- 
cal clusters of the present sample. The X-ray data are 
presented in form of contours of constant signal-to-noise 
above background in units of la, 2a, .... The optical data 
are from digitized Schmidt plates retrived from the 2nd 
version of the Digital Sky Survey (DSS2) archive. The 
clusters RS26, RS53 and RS75 in Fig.[Tll are not in the 
SDSS cluster catalogues of Gogo et al. (2002a) and Bahcall 
et al. (2003). 

All clusters of the present sample show a concentration 
of optical galaxies towards the X-ray center (in most cases 
clearly visible even in the less deep DSS2 images). The se- 
lection criteria described in Sect.Elthus appear quite use- 





Fig. 11. RS7 at 2 = 0.0581 with the X-ray flux S = 5.4 x 10"^^ ergs"^ cm''^ in the ROSAT energy range 0.1-2.4 keV. The 
cluster is not in the Bahcall et al (2003) sample. RS26 at z = 0.0663 with S = 5.1 x 10^^^ ergs~^ cm~^ is not in the Bahcall 
et al. (2003) and Goto et al. (2002a) catalogues. RS53 at 2 = 0.0425 with 5" = 3.7 x 10"^^ ergs"^ cm"^ is not in the Bahcall 
et al. and Goto et al. catalogues. RS75 at 2 = 0.1822 with S — 1.3 x 10~^^ ergs~^ cm~^ is not in the Bahcall et al. and Goto 
et al. catalogues. The cluster is also not in the RBF sample extracted from RASS-2 (Voges et al. 1999). 
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Fig. 12. Upper panel: Cumulative distribution of the neg- 
ative logarithmic Kolmogorov-Smirnov probability for a point 
source. The vertical line marks the threshold below which clus- 
ters are expected to be point-like in RASS-3. Lower panel: 
Cumulative distribution of the normalized differences between 
the observed hardness ratios to the theoretically expected hard- 
ness ratios for a 5keV cluster, normalized to the error of 
the hardness ratio measurement. The vertical line marks the 
threshold below cluster fluxes are expected to be significantly 
contaminated by AGN. All quantities are determined with the 
CCA method. 



ful: all X-ray cluster candidates passed the visual screening 
process. We have only excluded 7 targets which were ob- 
vious fragments of already detected objects. They could 
have been excluded already during the source detection in 
the maximum likelihood maps by introducing a minimum 
distance between two detected objects. However, from our 
past experience with NORAS and REFLEX we prefer to 
reject objects only after a careful visual screening at the 
end of the complete data processing when all data are in 
hand. 



7. Discussion and conclusions 

The present paper describes a method to search for clus- 
ters of galaxies by combining data obtained in different 
energy ranges. The main goal is to reach lower flux limits 
and to improve the reliability of a first identification of 
the sources. 

As an example, we combined X-ray data from RASS- 
3 and optical data from SDSS. The SDSS galaxy sample 
reaches much higher redshifts than RASS and can thus 
be used to guide a cluster detection in RASS down to 
lower X-ray flux limits. We could show that a detection 
of X-ray clusters in RASS-3 is possible down to a nomi- 
nal X-ray flux limit of 3 — 5 x 10~^'^ergs~^ cm~^ in the 
ROSAT energy band 0.1-2.4keV. This limit is about 5-10 
times deeper than REFLEX and gives a first qualitative 
impression how much deeper one can get with the new 
technique. 

In the following we discuss the results which are al- 
ready indicated by the analysis of the test area. Table 2 
shows that only 40 out of our 75 X-ray clusters (53%) are 
listed as X-ray sources in the ROSAT Bright and Faint 
Source (RBF) catalogues. The overlap can be increased 
possibly to 73% when a preliminary source list obtained 
with RASS-3 is used. It is thus not only the X-ray flux 
limit which can be decreased but also the X-ray complete- 
ness which can be increased when RASS and SDSS data 
are combined. 

At the nominal X-ray flux limit only 20% of the 75 
clusters show a signiflcant X-ray extent as measured with 
the GCA method (Fig.lEl Tab. 3). This smaU fraction 
is mainly caused by the comparatively small number of 
source counts (half of the sample has source counts be- 
tween 5 and 10). Usually, more than 30 source counts are 
necessary for the GCA method for a signiflcant detection 
of a source extent - the precise numbers do, however, also 
depend on the intrinsic extent of the source. 

One could worry about a signiflcant contamination by 
AGN. However, only 12% of the clusters have a known 
nearby AGN (Tab. 3) and only 5% of the X-ray clusters 
have an X-ray hardness ratio measured with GCA which 
is consistent with an AGN fFig. 1121) . The fraction of con- 
taminated clusters should thus be not very large, possibly 
around 10%. 

Before we discuss the direct comparisons with pub- 
lished SDSS-based cluster samples we want to note that 
it would not make much sense to test the number of 
SDSS clusters which are not present in our sample be- 
cause SDSS-based cluster surveys are expected to reach 
much higher redshifts and less rich systems compared to 
RASS-based cluster surveys. We are thus mainly looking 
for X-clusters which are not found in SDSS samples. 

The comparison of our X-ray cluster sample with pub- 
lic SDSS cluster samples yields the following results: 13 
of our 75 X-ray clusters (17%) are not in the catalogue of 
Goto et al. (2002a) and 23 clusters (31%) are not in the 
catalogue of Bahcall et al. (2003). At the comparatively 
bright X-ray fluxes S > 5.0-10"^'^ ergs~^ cm~^ where spu- 
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Fig. 13. X-ray luminosity vs. redshift for the total cluster sample (upper left), for clusters not in the Goto et al. (2002a) sample 
(lower left), for clusters not in the Bahcall et al. (2003) sample (upper right), and for the clusters neither in Goto et al. and 
Bahcall et al. (lower right). The redshifts used are estimated from the SDSS data with the matched filter described in Sect.|3| 



rious detections are unlikely, still 11 clusters (15% of our 
sample) are not in the Goto et al. sample and 19 clusters 
(25% of our sample) are not in the Bahcall et al. sample. 

However, when we combine the catalogues of Goto et 
al. and Bahcall et al. then the incompleteness of the opti- 
cal sample decreases to a level of 8%. The corresponding 6 
X-ray clusters which are not in the two SDSS-based cluster 
catalogues are marked as BG in the last column of Table 2 
and have the spectroscopic redshifts z — 0.0637, 0.0650, 
0.0663, 0.0425, 0.0206 and 0.1822. Most of these missing 
clusters (with one exception) are nearby and less rich sys- 
tems. This is also seen in the X-ray luminosity /redshift 
plots shown in Fig. ll3l In our X-ray cluster sample we have 
12 clusters with redshifts z < 0.07 and only 7 (58%) are 
in the combined optical catalogues of SDSS. One should 
note that the 6 X-ray clusters missing in the optical cata- 
logues on our test area of 140 square degrees could lead to 
an extrapolated number of 300 clusters over a total SDSS 
survey area of 7000 square degrees. We are thus discussing 
here a significant fraction of nearby systems. 

The difference between X-ray and optical cluster se- 
lection becomes apparent in Fig. 1141 where the redshift 
distribution of the total sample (unfilled histogram) and 
the distribution of clusters which are members of both 
the Goto et al. and the Bahcall et al. samples are com- 
pared. Ignoring the redshifts z > 0.3 where the sample 
of Bahcall et al. is incomplete (by construction), the red- 
shift distributions are significantly different as seen by the 
strong deficiency of nearby clusters in the combined Goto 
plus Bahcall catalogue. Whereas the cluster counts of the 
combined optical catalogue increase with redshift ~ which 
is typical for volume-limited samples - the X-ray sample 
shows a steady decrease with redshift - which is typical for 
flux-limited samples. Numerical simulations and compar- 



isons with deep XMM and Chandra observations will give 
further information about the final incompleteness level 
which can be reached with our method. 

We want to close the discussion of the sample with a 
few general remarks on the quite deep X-ray flux limit 
achieved from the combination of RASS and SDSS data. 
The resulting 'superfaint' flux limit - though optimal for 
the detection of high-redshift clusters and for compiling 
huge sample sizes - yields X-ray cluster candidates which 
consists in many cases of only 5-10 RASS X-ray source 
counts. The corresponding formal flux errors are of the 
order of 30 and 40%. A careful evaluation of the result- 
ing accuracies achievable in cosmological tests with such 
samples is thus necessary. 

In addition, more work is needed to better calibrate 
and test our method: (1) The analysis of a larger test area 
(Data Release 2) is expected to give significantly smaller 
statistical errors of the cluster number counts which is 
necessary for a better calibration of the likelihood thresh- 
old of the optical data. (2) Tests are under way which 
quantify the effect an optical sample like SDSS has on the 
selection of X-ray clusters and vice versa. In Popesso et 
al. (2004) the correlations between basic optical and X-ray 
cluster properties are studied and more tests will follow. 
One of their conclusions is that the z' SDSS photometric 
band and not the r' band used here gives the strongest 
correlations between optical and X-ray data, and should 
thus be the preferred band when we want to work with the 
smallest bias caused by including optical data in the detec- 
tion process. (3) The numerical simulations in Schuecker & 
Bohringer (1998) have to be performed under the specific 
conditions of the present survey in order to test selection 
effects caused especially by the small number of X-ray 
photons and surface brightness. (4) A direct comparison 
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Fig. 14. Redshift histogram of the clusters with spectroscopic 
redshifts. Filled steps show the redshift distribution of the clus- 
ters which are in common with Goto et al. (2002a) and Bahcall 
et al. (2003). 



with ROSAT, XMM-Newton, and Chandra pointed obser- 
vations will give model-independent estimates of the flux 
and completeness limits. 

7.1. Role of Virtual Observatories 

The present investigation gives a good illustration of the 
benefits achievable from the combination of high-quality 
optical and X-ray data. The combination of the two large 
databases allows us to fill the gap between nearby wide- 
angle and deep serendipitous X-ray cluster surveys with 
already existing data. The combination increases the com- 
pleteness of the samples and reduces the time-consuming 
spectroscopic follow-up observations. This method is also 
expected to be helpful when samples of the order of sev- 
eral lO'* X-ray cluster candidates expected from future 
wide-angle X-ray surveys like DUO or ROSITA must be 
identified. 

However, the quality and completeness of the cata- 
logues depends also on the model used to detect and char- 
acterize the galaxy clusters. In the present investigation 
only one set of values of model parameters (cluster core 
radius, cutoff radius, Schechter characteristic magnitude 
and faint end slope) was used to develope and illustrate 
the method. More realistic applications should test differ- 
ent sets of model parameters and use a finer sampling grid 
for cluster detection. In addition, cluster radii are known 
to increase with mass so that some kind of iterative pro- 
cedure is necessary to find the optimal model for each 
cluster. 

Further applications should also include the detection 
of clusters in so-called Sunyaev-Zel'dovich maps obtained 
in the centimeter or millimeter range. Such kinds of maps 
will be provided by, e.g., the PLANCK mission, and will 
surely increase the reliability of the resulting catalogues. 



We are thus left with the conclusion that the quality, 
the completeness limit, the detection limit, and the tele- 
scope time for source identification can be optimized when 
enough computer power to sample the high-dimensional 
parameter spaces is available. Moreover, a fast and easy 
handling of data is needed to work with very large 
databases probably located at different physical nodes. 
Virtual observatories like the GAVO provide the infras- 
tructure for such kind of investigation and are surely nec- 
essary to make full and efficient use of already existing 
and future observed data. 

GAVO provides an all-sky X-ray event file including 
many characteristics of the invididual photons. The pho- 
tons are retrieved from the archive with a general cone- 
search algorithm. We are currently in the process to imple- 
ment the cluster search described here within the GAVO 
grid. This final step is crucial for the analyses of larger sur- 
vey areas with more cluster models and a finer sampling 
grid of the detection cells. This will allow us to analyse in 
a much more efficient manner the SDSS Data Release 2 
sample which will become available soon and RASS-3. 

The method described in the present investigation 
needs about 1 hour computing time per square degree 
on a standard workstation with one processor. Running 
the method in a controlled manner over a grid of many 
computers is thus essential for a proper analysis within 
reasonable times. Our plans include the provision of the 
cluster (source) detection algorithm as a service of GAVO 
so that the procedures can be used in a user-friendly man- 
ner by a larger community of interested researchers. 
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Tab. 2. Catalogue of RASS-3/SDSS cluster candidates. Col. 1: Cluster number. Col. 2-3: cluster coordinates in degrees. Col. 4: 
column density of neutral galactic hydrogen in units of 10^" cm^'^. Col. 5: unabsorbed X-ray flux in units of lO^^"* ergs~^ cm~^ 
in the ROSAT energy band 0. 1-2.4 keV. Col. 6-7: detection likelihoods on RASS-3 photon maps and SDSS galgixy maps. Col. 8- 
9: estimated cluster redshifts from RASS-3 and SDSS data. Col. 10: Spectroscopic cluster redshift. Col. 11: Number of cluster 
galaxies (after rejection of outliners) used for the determination of the spectroscopic cluster redshift. Col. 12: Note whether the 
cluster is not found in the RASS source catalogues (V), and not found in the SDSS cluster catalogues compiled by Bahcall et 
al. (B, 2003) and Goto et al. (G, 2002a). For more details see main text. 
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-0.7342 


7.88 


9.35 


12.491 


51.521 


0.34 


0.18 








G 


RS73 


53.9795 


-0.6320 


8.01 


5.58 


4.455 


20.131 


0.22 


0.22 


0, 


.2741 


2 


V23 


RS74 


55.6623 


-0.2708 


7.79 


6.83 


5.461 


42.824 


0.28 


0.40 








V2B 


RS75 


56.3280 


1.0026 


10.60 


12.86 


6.229 


27.567 


0.08 


0.16 


0, 


.1822 


4 


V2BG 
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Tab. 3. Catalogue of RASS-3/SDSS cluster candidates. Col. 1: cluster number. Cols. 2-3: cluster coordinates in degrees. Col. 4: 
source counts. Col. 5: countratcs from matched filter and formal la errors. Col. 6: countrate from Growth Curve Analysis GCA 
and formal la errors. Col. 7: normalized hardness ratio (see text). Col. 8: negative logarithmic Kolmogorov-Smirnov probability 
for a point source. Col. 9: optical richness. Col. 10: further remarks (CONT: AGN within 1 arcmin from cluster position in NED; 
CUT: cluster only partially in test area of matched filter). 



No. 


R.A. 


DEC. 




A 


ctr(GCA) 

J 


HR 


EXT 


Riclmess 


Remarks 




2000.0 


2000.0 




fs-^l 

L*^ J 


GCA 


GCA 


SDSS 




RSOl 


0.3581 


-0.0312 


12.5 


n'^fi+n mo 


yj ^yj ^ A i_ \j A-Kj 


0.59 


1.4 


79.0 




RS02 


4.3972 


-0 8836 


11.3 


033±0 010 


0.023+0.011 


0.75 


0.4 


90.3 


ZwCl 0015.1-0108 


RS03 


5.0566 


0.1015 


13.9 


0.036±0.010 


0.026+0.009 


-1.46 


0.2 


128.0 


RXC00201.1 


RS04 


5.1680 


-0 9858 


8.8 


025±0 008 


018+0 008 


-1.52 


0.0 


41.4 


2nd CLST z = Q 1071 

^l.X\X l—iV-J X \J m A^\J 1 X 


RS05 


5.1943 


-1.0844 


6.0 


017±0 007 


015+0 007 


0.57 


2.1 


37.0 




RS06 


5.4073 


-0.8459 


8.0 


021±0 007 


015+0 007 


-0.27 


1.4 


113.3 


A0023 


RS07 


5.6274 


-0.9570 


10.2 


026±0 008 


026+0 015 


0.80 


0.0 


52.3 


zn =Q 12921 


RS08 


5.7515 


-0.1543 


12.1 


030±0 009 


022+0 010 


0.02 


0.1 


151.3 


A0025 


RS09 


5.7988 


0.1977 


14.5 


0.038±0.010 


0.022+0.010 


-0.30 


1.0 


62.1 


ZG = 0.25402 


RSIO 


7.1309 


-0.0857 


9.1 


019±0 006 


013+0 006 


0.51 


0.6 


46.8 


ZwCl 0025.9-0.021, CONT 


RSll 


7.3730 


-0 2509 


39.7 


081-1-0 013 


082+0 018 


-1.11 


1.5 


92.0 


ZwCl 0027 0-0036 

£-J w v_'i yjyjij 1 .yj yjyjtjyj 


RS12 


8.0414 


-0.6574 


9.7 


0.015±0.005 


0.009+0.005 


-0.54 


0.2 


56.5 


ZG = 0.24267 


RS13 


8.5759 


0.8777 


13.6 


032±0 009 


0.021+0.011 


-0.95 


1.0 


94.2 




RS14 


9.3587 


0.1414 


5.8 


019±0 008 


017+0 009 

\j mKj A- i _i_ yj *yjyj 


-4.74 


0.3 


54.0 


CONT 


RS15 


10 7036 


0.2149 


12.2 


039±0 Oil 


0.027+0.012 


0.57 


0.2 


122.2 




RS16 


10 9596 


0.1187 


8.1 


026±0 009 


026+0 015 


-2.01 


2.7 


74.3 


ZwCl 0041 1-0008 


RS17 


11.4570 


-0.8486 


21.3 


067±0 015 


068+0 017 


-0.38 


1.1 


51.8 


A0095 


RS18 


11.5751 


0.0167 


10.0 


0.033±0.010 


0.028+0.011 


-0.85 


0.1 


83.3 




RS19 


11.6170 


-0.2219 


9.5 


031±0 010 


0.028+0.011 


0.16 


0.3 


62.9 


zn — 26537 


RS20 


1 1 8959 


-0.8947 


10.5 


035+0 01 1 

\J . \J fJ ' J 1 \J • \J A. A. 


020+0 009 

yj • yjiuyj i yj ,yjyj u 


0.28 


1.0 


49.6 


AOlOl 


RS21 


13.9519 


0.6488 


29.1 


0.078±0.014 


0.059+0.015 


-0.30 


2.7 


67.1 


A0116 


RS22 


14.0756 


-1.1735 


207.3 


621±0 043 


1 413+0 071 


0.41 


-9.9 


56.1 


A0119 CUT 


RS23 


15.6415 


-0.8327 


7.3 


019±0 007 


015+0 007 

• vy X ty — i— yj *yjyj 1 


0.14 


0.0 


74.6 




RS24 


15 6763 


1.1263 


31.9 


077±0 014 


060+0 014 

yj • yj\jyj -A^\j *yj 


0.50 


1.4 


128.8 


CONT 


RS25 


15 9883 


-0.8067 


6.6 


016±0 006 


015+0 015 


-0.50 


1.5 


64.7 




RS26 


16.0781 


-0.7695 


10.3 


025±0 008 

y-} • \J^^ ^^\J »yjyjKJ 


012+0 006 


-2.98 


0.5 


47.5 


RASSCALS 087 


RS27 


16.0791 


-0.4288 


16.8 


0.040±0.010 


0.023+0.008 


-0.02 


1.1 


101.6 


CONT 


RS28 


16.2231 


0.0484 


26.4 


0.064±0.012 


062+0 013 

yj M yj\j —1— yj *yj A-f^ 


-2.45 


2.5 


174.5 


ZwCl 0102.4-0012 


RS29 


1 6 3765 


-0.1699 


14.2 


035-1-0 009 


030+0 019 


-0.60 


0.9 


64.0 


2ri — 231 33 


RS30 


16.5541 


0.8009 


4.1 


0.122±0.019 


0.012+0.006 


0.11 


2.4 


90.0 




RS31 


16 7056 


1.0513 


69.3 


165±0 020 


0.149+0.020 


-0.55 


0.4 


84.1 


RXCJ 0106 8 CONT 


RS32 


16 7805 


-0 3280 


13.4 


032±0 009 


032+0 015 

yj M yj\j ^ —i-^ yj ^yj A.^j 


0.12 


0.8 


139.2 


zn — Q 26537 


RS33 


17.8626 


0.5741 


6.9 


0.017±0.006 


0.015+0.007 


0.37 


0.3 


104.2 


ZG = 0.36748, CONT 


RS34 


18.7749 


3502 


188.7 


500±0 036 


435+0 039 


-3.35 


24.9 


69.6 


A0168, RXCJ0114.9 


RS35 


20.5041 


3503 


24.5 


049±0 010 


068+0 020 

\j *\j\jyj .^^yj m\j 


0.82 


2.7 


86.9 


A0181, RXCJ0121.9 


RS36 


21.8232 


0.3523 


27.3 


0.065±0.012 


0.059+0.013 


0.29 


2.4 


181.7 


ZG = 0.35613 


RS37 


22 0578 


-0.6777 


8.4 


020±0 007 


0.023+0.011 


-2.47 


0.2 


66.5 


zn — 40152 


RS38 


22 6058 


0.4817 


6.1 


014+0 006 

yj • yj 1 1^. yjyjyj 


01 1 +0 006 

yj .yj A. A — L \j .yjyjyj 


-1.39 


0.6 


107.2 




RS39 


22.9211 


0.5533 


39.1 


0.095±0.015 


0.090+0.018 


-1.87 


5.9 


53.1 


A0208, RXC0131 


RS40 


23.7341 


-0.9379 


6.6 


0.016±0.006 


0.016+0.007 


-1.70 


0.2 


46.0 




RS41 


24.3671 


-0.1861 


7.0 


0.023+0.009 


0.021+0.010 


0.62 


0.0 


162.1 


ZG = 0.36748 


RS42 


26.6947 


-0.6782 


39.3 


0.099+0.016 


0.081+0.015 


-4.09 


0.1 


34.6 


CONT 


RS43 


27.2777 


-0.4271 


19.6 


0.049+0.011 


0.027+0.015 


0.29 


0.5 


64.0 


ZG = 0.33344 


RS44 


27.3363 


-1.1850 


6.7 


0.016±0.006 


0.023+0.010 


-0.70 


0.1 


160.9 


ZG = 0.29940 


RS45 


27.3458 


-0.3773 


15.3 


0.038±0.010 


0.011+0.006 


0.16 


0.2 


81.7 


ZwCl 0146.8-0035 


RS46 


27.8459 


-1.0133 


6.2 


0.015+0.006 


0.012+0.006 


-0.94 


0.2 


44.3 




RS47 


28.1781 


1.0044 


85.8 


0.206+0.022 


0.209+0.027 


-1.93 


3.6 


180.5 


A0267, RXCJ0152.7 


RS48 


28.6185 


-0.6531 


9.7 


0.023±0.007 


0.017+0.007 


0.14 


0.0 


108.9 


ZG = 0.43556 


RS49 


29.0463 


1.0481 


23.3 


0.058+0.012 


0.065+0.016 


0.14 


4.2 


96.2 


A0279 
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